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O ' Abstract 
^ I 

■ The reasons for using low-lying Dirac eigenmodes to probe the local structure of topo- 

^ , logical charge fluctuations in QCD are discussed, and it is pointed out that the qual- 

I itative double-peaked behavior of the local chiral orientation probability distribution 

^ ■ in these modes is necessary, but not sufficient, for dominance of instanton- like fluctua- 

^ . tions. The results with overlap Dirac operator in Wilson gauge backgrounds at lattice 

I spacings ranging from a ~ 0.04 fm to a ~ 0.12 fm are reported, and it is found that 

' the size and density of local structures responsible for double-peaking of the distribu- 

^> . tion are in disagreement with the assumptions of the Instanton Liquid Model. More 

^ , generally, our results suggest that vacuum fluctuations of topological charge are not 

' effectively dominated by locally quantized (integer-valued) lumps in QCD. 

1 Introduction 



It has been recently suggested [1] that questions about the dynamical nature of topological 
charge fluctuations in the QCD vacuum are worth revisiting. The notion that gauge field 
topology is relevant for understanding QCD started with the discovery of instantons [2], and 
their subsequent use as a basis for, among other things, the qualitative resolution of the Ua{^) 
problem [3], and the discovery of 9-dependence of QCD physics [4]. These successes were as- 
sociated with the use of semiclassical methods and a concrete /oca/ picture of the vacuum, the 
instanton gas picture [5], characterized by formation of well-separated (anti)self-dual lumps 
of quantized topological charge. However, it was soon realized that such a vacuum does 
not lead to confinement, thus putting the relevance of instantons immediately into question. 
Refiecting upon this situation, Witten has argued [6] that large quantum fiuctuations asso- 
ciated with the confining vacuum naturally generate large fiuctuations of topological charge. 
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and lead to similar qualitative effects as those usually ascribed to instantons. In fact, he has 
conjectured that instantons are not important dynamically, arguing that the semiclassical 
picture is invalidated by large quantum corrections, and suggesting that topological charge 
fluctuates in a more or less continuous manner. 

Nevertheless, QCD instantons remained rather popular as a frequently preferred (and 
analytically accessible) way of thinking about vacuum topology. The instanton solution was 
used as a basis for developing a rather successful semiclassically-motivated phenomenology, 
the Instanton Liquid Model (ILM) [7], where correlations among instantons were introduced 
to suppress the infrared divergences present in the instanton gas. Even though termed a 
"liquid" , the corresponding vacuum is quite dilute with (anti)instantons of radius p ~ 5 fm 
and density n ~ 1 fm~^ occupying a small fraction of space-time volume and retaining 
their identity. This setup allows for an interesting mechanism of spontaneous chiral sym- 
metry breaking [8], in which the t'Hooft "would-be" zeromodes associated with individual 
(anti)instantons mix, and supply the finite density of near-zeromodes required by the Banks- 
Casher relation [9]. While the ILM is a phenomenological model, this elegant mixing picture 
invokes the impression that instantons play an important microscopic dynamical role in the 
QCD vacuum. A primary aim of our investigation is to examine whether the microscopic 
relevance of ILM picture can be justified. 

The support for the foundations of the ILM has frequently been drawn from lattice 
QCD simulations, using equilibrated gauge configurations as typical representatives in the 
path integral (see e.g. Ref. [10]). Indeed, it is quite reasonable to expect that lattice QCD 
will eventually provide us with detailed answers about the nature of topological charge 
fiuctuations. However, finding a clean and satisfactory way to infer this information from 
lattice QCD has proven to be a nontrivial issue due to the fact that lattice gauge fields 
are found to be very rough at the scale of the lattice spacing. The apparent necessity to 
eliminate the short-distance fluctuations in some way resulted in manipulating the gauge 
fields in various cooling or smoothing procedures. While useful for estimating the global 
topological charge, these techniques are significantly biased as to the local structure of the 
vacuum they reveal, and there is always an inherent subjective element present. 

The idea of using fermions to study global topology is relatively old [11], but was not 
very much exploited due to both computational demands and problems with chiral symmetry. 
However, the suggestion to use low-lying fermionic modes to extract the information about 
/oca/ fiuctuations of topological charge is recent [1] (see also Ref. [12]). In the light of the 
above remarks, the preference for using fermions in this respect is quite reasonable, but 
there is also a good physical motivation for doing so. In particular, local fluctuations of 
topological charge are of interest because they are believed to cause hght quarks to generate 
the rj' mass [3, 6] and the chiral condensate [8], for example. It is thus physically very 
natural to look for the imprint of these fiuctuations in the low-lying Dirac eigenmodes which 
dominate the propagation of light quarks. 

The concrete proposal of Ref. [1] is to investigate the behavior of local chirality in the low- 
lying modes. This approach uses the fact that the instanton-like gauge fluctuations would 
leave a specific imprint in the individual modes, thus offering an opportunity to examine the 
consistency of the instanton picture. More precisely, it was argued that if there are extended 
regions of (anti)sclf-duality in the gauge background, then the probability distribution of 
the chiral orientation parameter X (X-distribution) over the regions of strong field should 
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exhibit double-peaked behavior in the vicinity of extremal chiralities. The initial study of the 

X-distribiition in Rcf. [1] indicated a flat behavior and thus a qualitative discrepancy with 
the instanton picture. However, later studies on finer lattices and/or with chirally symmetric 
fermionic actions [13, 14] revealed a visible degree of double-peaking. This situation precludes 
any definite conclusions on purely qualitative grounds since the double-peaked structure is 
necessary but not sufficient for instanton dominance. It is a main purpose of this work to 
offer a more quantitative point of view. 

In the first part of the paper (Sec. 2) we explain in detail why we choose to rely on the in- 
formation encoded in low-lying Dirac modes, rather than using the conventional approaches 
for probing the gauge field topology on the lattice. We then concentrate on clarifying the 
possible implications of the X-distribution and argue that a double-peaked structure is not 
uniquely associated with self-duality. For example, if the gauge fields form approximately 
quantized ^ isolated lumps of topological charge, then strong double-peaking is expected as 
well, even when such lumps are not self-dual. Also, some chiral peaking in the low eigen- 
modes may occur even for non-quantized topological charge excitations through some as-yet- 
unknown mechanism. The observed behavior of the X-distribution is thus not necessarily 
associated with an instanton-dominated vacuum. We suggest that the tentative conclusions 
of such nature be supplemented by quantitative comparison of QCD X-distributions to those 
for the existing ILM ensembles. The material in Sec. 2 extends and complements the general 
discussion of Ref . [1]. 

In the second part of this paper (Sec. 3), we present new results for X-distributions from 
eigenmodes of the overlap operator in Wilson gauge backgrounds at lattice spacings ranging 
from a ^ 0.04 fm to a ^ 0.12 fm. In an attempt to verify whether the observed behavior 
can be ascribed to gauge structures with the parameters of the ILM, we determine the size 
and density of fermionic structures that contribute to the peaks of the X-distribution. This 
is motivated by the fact that for instanton-like excitations there is a direct correspondence 
between gauge and fermionic structure. Our consistency check leads to significant disagree- 
ment with the ILM. In fact, the lattice spacing dependence of the average radius (defined 
without reference to the t'Hooft zero mode profile which does not fit reasonably) indicates a 
finite continuum limit with value i? O.lSfm, while the density of structures at a 0.04fm 
is n ~ 50 fm~^, and we cannot exclude the possibility that it diverges in the continuum limit. 
We then point out that this not only disagrees with the ILM, but also suggests that the bulk 
of topological charge in QCD is not locally quantized in approximately integer units. We 
give several arguments supporting this interpretation. Needless to say, if confirmed, this 
would have profound implications for the possible microscopic explanation of spontaneous 
chiral symmetry breaking. Finally, we give several arguments that our data is not contam- 
inated by lattice artifacts usually referred to as "dislocations", and support this conclusion 
by presenting the data for the Iwasaki gauge action. A preliminary version of this work is 
presented in compact form in Ref. [15]. 

^When we refer to local quantization of topological charge we always imply quantization in units of ±1. 
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2 Cooling or Smoothing vs Fermions 



The structure of equilibrium lattice Monte Carlo configurations provides a unique window 
for examining the nature of gauge fluctuations in the QCD vacuum. It also represents a 
distinctive way of thinking about how low-energy phenomena in QCD arise in terms of 
fundamental degrees of freedom. For example, as wc will examine in some detail below, in 
the ILM mechanism of spontaneous chiral symmetry breaking the condensate arises due to 
mixing of t'Hooft's "would be" zeromodes associated with (anti)self-dual lumps, carrying 
approximately quantized topological charge. This line of reasoning inherently assumes that 
this happens at the "configuration level" with the vacuum dynamically generating the gauge 
potentials with these properties. If this is the case then the picture has a fundamental 
microscopic meaning. Otherwise, it represents only a phenomenological modeling. 

Unfortunately, examining lattice gauge fields directly leads to ambiguous results, mainly 
due to the fact that fields are typically rough even at the scale of a single lattice spacing. 
While the situation will improve with the approach to the continuum limit, this is probably 
not so when comparing the behavior at fixed physical distance. It is thus perhaps inevitable 
that some sort of filtering procedure be used to interpret the local structure of lattice gauge 
fields. In this section, we examine various ways of approaching this issue. 

Before proceeding to discuss this in more detail, it is useful to fix the language that will be 
used in what follows. In particular, we will use the term "isolated lump" or simply "lump" to 
denote a local structure in the gauge field that can be enclosed by a hypersurface on which the 
field is approximately pure gauge. Consequently, the lump contains approximately integer- 
valued topological charge. A tendency towards local "lumpy structure" is one of the inherent 
properties of the ILM picture which distinguishes it from the situation when the gauge field 
fluctuates very inhomogeneously, forming peaks of topological charge density, but without 
any tendency for such quantization. 

2.1 Cooling 

One popular way of addressing the problem of short-distance fluctuations is the cooling 

method [16]. We will be rather generic when using this term, assuming just that it is a 
local minimization procedure for the gauge action with initial state being the Monte-Carlo 
generated QCD configuration. We will not be concerned with various implementations or 
variations on the original idea. 

With the elementary step being local, it is naturally expected that if the original config- 
uration can be assigned a global topology in some way, then it will be preserved to a large 
extent during such a procedure. One can thus turn the argument around and to define the 
global topological charge of the original configuration through field-theoretic definition, let's 
say, on the corresponding cooled counterpart. Such a definition may not be entirely satisfac- 
tory (the number of cooling sweeps is rather subjective and non-unique) but it can be made 
into a fairly well-motivated working scheme especially for the improved versions [17]. 

However, the situation is different if one is interested in the local structure of topolog- 
ical charge fluctuations. In the continuum, there are true local minima of the action for 
strictly self-dual or strictly anti-self-dual fields. The superpositions of not overwhelmingly 
overlapping instantons and antiinstantons, while not exact minima, correspond to plateaus 
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(or shallow valleys) in the action profile. After a few cooling sweeps when the gauge field 
undergoes large changes, the local minimization can bring the configuration to the vicinity 
of such a plateau. Here the fields become naturally smooth, the field-thcorctic definition of 
topological charge leads to values that cluster around integers, and the evolution in the con- 
figuration space slows down. Further cooling possibly leads the configuration into another 
plateau, and eventually into the global minimum with nonperturbative fields completely re- 
moved ^. By the nature of the argument, in the vicinity of the plateau the configuration will 
necessarily resemble a (multi)instanton-antiinstanton state regardless of the local properties 
of the original configuration. In other words, a possible observation of lumpy structure and 
local (anti)self-duality in the cooled configuration is to be expected, and cannot be used as 
an independent, logically satisfactory input for making conclusions about local fiuctuations 
of topological charge or about the dynamical importance of instanton-like gauge fiuctuations 
in the QCD vacuum. 

To illustrate this point, it is perhaps instructive to imagine that we are interested in 
studying the local structure of a sufficiently strongly-coupled lattice QCD vacuum, where 
instanton-like fluctuations are not expected to play any role. Consequently, an unbiased 
approach should find no traces of them. However, cooling the corresponding equilibrium 
configurations will still lead to plateaus and locally lumpy, self-dual behavior. Can we use 
this as a basis to conclude that instantons play a significant role in the strongly-coupled 
lattice QCD vacuum? Certainly not. 

2.2 Smoothing 

The inherent reason why cooling is a biased way of studying local properties of gauge fiuctu- 
ations is that it uses the gauge action as a basis for the procedure. This is not easily curable 
by improvements. However, one can also attempt to eliminate the unphysical short-distance 
fiuctuations by some sort of smoothing without reference to the action. Cooling itself is a 
smoothing procedure in the sense that it replaces the original configuration with a smoother 
one. For the discussion in this section, however, we will be more specific in that respect, and 
call smoothing any procedure that can be interpreted as an action-independent space-time 
averaging of the fields. The prototype for this would be for example an APE-smearing [18] 
that was actually used in this context [19]. 

While smoothing is a less invasive approach than cooling, the subjective element in 
determining just how much smoothing is enough remains. This is troublesome because 
smoothing, by its nature, is also biased towards lumpy structure and can qualitatively change 
the local behavior. To see that, consider a very rough, inhomogeneous gauge field with 
various structures in it. Upon smoothing, the field varying over short distances will be quickly 
removed or reduced to a smooth four-dimensional "bump" if there is an underlying long 
distance structure present. Also, the lower-dimensional structures, such as ridges, sheets etc, 
will be eliminated soon by four-dimensional averaging. Thus, as the smoothing progresses, 
one naturally expects the stage when the configuration will be considerably lumpy with 
pliysieal fields eoiK eiitratofl in tlic lumps. Further smoothing will cause the lumps to grow 

^For improved cooling it is in principle possible that the configuration remains indefinitely in the vicinity 
of local self-dual minimum (with instanton-antiinstanton pairs annihilated). 
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in size, to overlap, and eventually to reach a homogeneous situation with all physical fields 
removed. An unsettling question is where in this process one should stop and claim that the 
fields at the given stage represent the filtered local structure of the original configuration. 

2.3 Fermions 

Fortunately, there is a clean and meaningful strategy for approaching the above issues, 
namely to rely on the fermionic response to the corresponding gauge background. To explain 
that, it is easiest to first think in terms of a theory in the continuum where a typical 
configuration is also expected to have a lot of ultraviolet fiuctuations. Using fermions is 
quite plausible since we know that, for sufficiently smooth gauge fields, fermions refiect 
global topology exactly [20] and for non-diffcrcntiabic gauge fields the index of the Dirac 
operator can actually serve as an extended definition of topological charge. It is thus natural 
to expect that the local structure can be inferred from fermionic response as well, e.g. by 
studying the divergence of the fiavor-singlet axial- vector current. 

The new idea which makes the fermionic approach attractive is to look for the imprints 
of topological charge fluctuations in the low-lying eigenmodes of the Dirac operator [1]. This 
is physically well-motivated, practical and, at the same time, has the potential to naturally 
solve the problem of ultraviolet fluctuations without the subjective element involved. Indeed, 
the space-time structure of the low-lying modes is naturally smoother than that of the gauge 
fields themselves. This fact has two origins. First, the analogous Schrodinger-hke eigenvalue 
problems typically yield smooth stationary states even when the underlying potential has 
discontinuous jumps. Secondly, among the eigenstates, the infrared modes are expected to 
be least sensitive to the short- distance features on top of the long-distance structure in the 
potential^. At the same time, it is necessary to realize that the underlying physics we want to 
understand is the mechanism of how gauge fiuctuations cause the light quarks to propagate 
in such a way as to form a quark condensate and to give the r]' its mass. In other words, the 
goal is to understand the dynamics underlying the propagation of light quarks. This physics 
is encoded in the low-lying Dirac modes [1], and this is why we say that concentrating on 
these modes ( "fermion filtering" ) is physically well motivated. 

One possible approach is to infer the information about topological charge fiuctuations 
from the individual low-lying modes [1]. This method is indirect and is motivated by the 
instanton picture of the vacuum. In other words, the properties of instant on- like gauge 
fiuctuations imply a specific qualitative behavior of the individual modes, thus giving an 
opportunity to check whether the observed structure is consistent with that scenario. The 
specific proposal of Ref. [1] is to study the probability distribution of the local chiral orien- 
tation parameter X{n) defined by 



over the points n where the cigcnmodc (ip^ipn) is large. Here 2-?/^^ = 2iI}r = (1+75)-?/' 

are the left and right components of the eigenmode ip. X{n) is a local angle in the \ijjL\- 
IV'rI plane rescaled so that X{n) ~ —1 for purely right-handed and X{n) — +1 for purely 

^Low-lying modes remain sensitive to the isolated small structures as they should. 




(1) 
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left-handed spinor For ipL,fpR generated in a random independent fashion, the local 
orientation parameter would be uniformly distributed between —1 and +1. 

2.4 The Implications of the X-distribution 

The usefulness of the X-distribution was discussed in detail in Ref. [1], where it was argued 
that if the gauge background contains extended regions of strong (anti) self-dual fields, then 
the X-distribution should exhibit peaks near the extremal values. This is a consequence of 
the fact that in the eigenvalue problem for the Dirac operator, the (anti)self-dual part of the 
gauge field enters as a potential term for the (right) left component of the eigenmode. 

Here we would like to emphasize that this is not a unique way for the double-peaked 
behavior to arise. To illustrate that, consider subjecting the fermion to a background con- 
figuration that is lumpy in the sense that we have defined it here, i.e. consisting of relatively 
isolated lumps surrounded by approximately field-free regions. (Such configurations can be 
artificially prepared irrespective of whether they resemble typical equilibrium configurations 
of QCD.) The fields comprising the lumps have no definite duality properties but yet, carry 
approximately quantized topological charge. A particular lump can thus be thought of as 
an individual object and, in the absence of all the other lumps, it would induce zero modes 
of the Dirac operator with chirality dictated by the index theorem. For simplicity, consider 
a situation with just two lumps with topological charges Q — ±1. The residual interaction 
will cause these "would-be" zero modes to mix and, to a first approximation, the true eigen- 
modes would be the two linear combinations of the strictly chiral "would-be" zeromodes 
localized on the lumps. In the resulting topological near-zero modes, the X-distribution will 
be strongly peaked at ±1. 

What we have described above is the same scenario that is a basis for the ILM mechanism 
of spontaneous chiral symmetry breaking [8] . The point is that self-duality is not essential for 
this argument since we didn't need to invoke instantons at all. It is only the local quantization 
of topological charge that matters and yet, the peaked X-distribution is expected, as well as 
spontaneous chiral symmetry breaking if such configurations turned out to be dynamically 
important. 

We emphasize this point to show that the qualitative double-peaked behavior of the 
X-distribution alone does not in itself provide a basis for the verification of the instanton 
picture. It is a necessary but certainly not a sufficient condition for such a conclusion. Apart 
from the above argument, there can possibly be other mechanisms for producing the peaked 
distribution (as well as spontaneous chiral symmetry breaking). In fact, we will argue later 
that there are reasons to believe that topological charge is actually not locally quantized in 
QCD. Nevertheless, the study of quantitative characteristics of the X-distribution can still 
be very useful for distinguishing various scenarios as well as for other purposes [14] . 

2.5 Lattice Fermions 

When considering the lattice-regularized theory we have the usual extra freedom in choosing 
the lattice action. According to the standard universality assumption, all local actions with 
appropriate symmetries and correct classical behavior should give consistent results suffi- 
ciently close to the continuum limit. However, comparing results from different discretiza- 
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tions at finite lattice spacing usually requires some care. When asking questions about the 
nature of fluctuations in the pure gauge vacuum, it is natural to fix a particular lattice 
gauge action and study the theory at different lattice spacings to be able to extrapolate to 
the continuum limit. Moreover, if we choose to examine the nature of these fluctuations 
using a lattice fermion, then it is essential that the fermionic action be kept flxed as well, 
even though it only plays the role of a probe. 

To illustrate this point, consider some lattice Dirac operator D{U) and the set of related 
operators Dn{U) = D{Un), where f/„ is the gauge configuration obtained from U by per- 
forming n smearing steps (e.g. APE steps) with other parameters of the smearing procedure 
flxed ^. If D{U) is an acceptable lattice Dirac operator, then Dn{U) is acceptable as well, 
assuming that C/„ exists. However, if we investigate the local properties in the equilibrium 
ensemble {U}eq using then the answers we get will clearly depend on n. In particular, 
-Diooo will probably indicate much larger structures in the gauge field than Di. Neverthe- 
less, one is in principle allowed to use i^iooo as long as it is used consistently in meaningful 
extrapolations to the continuum limit, which might be difficult. However, one can hardly 
conclude anything meaningful about the behavior in the continuum limit by comparing the 
results from Di at one lattice spacing and from i^iooo a different lattice spacing. 

Regarding the choice of the fermionic action, we would like to point out one particular 
aspect illustrated by comparing the two extremes represented by the Wilson-Dirac operator 
and the overlap operator based on it. The overlap operator has exact lattice chiral symmetry 
at finite lattice spacing, allowing for continuum-like theoretical analysis with respect to chiral 
symmetry and topology. On the other hand, it is non-ultralocal and, due to extended gauge 
connections, effectively samples the gauge field potential over nonzero physical distance ^. 
This should not matter for the continuum limit, based on the expectation that the operator is 
local over equilibrium ensembles of Wilson pure gauge theory at sufficiently weak coupling. 
The Wilson-Dirac operator, on the other hand, has mutilated chiral properties and when 
using it in this context, one relies heavily on the assumption that this problem will go 
away in the continuum limit. At the same time, the Wilson-Dirac operator has a perfect 
resolution in the sense that there is no additional smoothness beyond the fact that we end 
up inspecting the infrared eigenmodes. It would be very interesting to know whether the two 
actions provide for consistent continuum extrapolations for quantities related to topology. 
In this work, we use an overlap operator which, even though expensive to implement, is very 
convenient for its theoretical advantages. 

The above considerations also suggest that there may be a more general complementar- 
ity between the degree of chiral symmetry and the resolution for lattice fermionic actions. 
Indeed, if we start with a maximally ultralocal operator, we can try to improve its chiral 
properties by including couplings at larger distances and adding more complicated gauge con- 
nections, but by doing so we worsen the resolution. On the other hand, if one starts with the 
operator with exact lattice chiral symmetry, then this operator must be non-ultralocal [22], 
and an attempt to improve its resolution by dropping couplings at large distances will result 
in deteriorated chiral properties. 

^See Ref. [12] for example. 

''Due to this "chiral smoothing" the features smaUer than some physical threshold will not be resolved at 
finite lattice spacing. The results of Ref. [21] are a manifestation of this fact. 
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3 Low-Lying Modes of the Overlap Operator 



We have calculated and examined the low-lying modes of the overlap operator in Wilson 
gauge backgrounds over a wide range of lattice spacings. The massless overlap operator [23] 

D^p l + {Dw-p)-r = (2) 

L ^{Dw-p) + {Dw-p)\ 

with D\Y the Wilson-Dirac operator and p = 1.368 = 0.19) is used throughout this paper. 
The rational approximation for the matrix sign function [24, 25] was utilized to implement D 
and a small number of eigenmodes of j^D^ were projected out to both increase the accuracy 
of the approximation and to speed up the convergence. 

We used the Ritz variational method [26] to obtain the low-lying eigenmodes of = 
D^D. As a consequence of normality {[D,D^] = 0) and 75-hermiticity {D^ = 75/^75 ), 
the nonzero low-lying modes of are doubly degenerate and the eigenmodes of D can be 
constructed in the corresponding subspaces since [H^, D] = 0. Also, is proportional to the 
chirally non-symmetric part of D, implying that [75, H'^] = 0. This allows for diagonalization 
in separate chiral sectors in different runs, thus speeding up the process and easing the 
memory requirements for large lattices. The typical accuracy of calculated eigenvalues (as 
measured by differences of complex-conjugate pairs) is about one part in 10^ or better. 

To avoid confusion, it is useful to fix the language that will be used when discussing the 
local properties of a given eigenmode -0 of D. We will refer to 

d{n) = i^ni^n c(n) = iJnl5'ipn (3) 

as density and chirality respectively, and to X{n) of Eq. (1) as chiral angle. Since d{n), 
c{n) and X{n) are identical for eigenmodes corresponding to complex conjugate nonzero 
eigenvalues — Jbipx), we will treat ipx, as a pair. 

The parameters of the Wilson gauge ensembles used are hsted in Table 1. Configurations 
for the three finest lattices arc separated by twenty thousand sweeps. The quoted values of 
lattice spacings were obtained from the Sommer parameter using the interpolation formula 
given in Ref. [27]. Our finest lattice spacing is outside the interpolation range and we linearly 
extrapolate from P — 6.5. Note that the physical volumes involved are such as to contain 
on average 3-4 (anti)instantons if the ILM scenario is relevant, ensuring that the mixing of 
't Hooft zero modes would take place. For all the configurations we have calculated the 
eigenfunctions for the zeromodes and at least two pairs of near-zeromodes. 
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Table 1: Ensembles of Wilson gauge configurations. 
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Figure 2: Dependence of < |X| >, Eq. (4), on the lattice spacing. The quoted uncertainty 
represents a rough estimate using the asymmetry \pxi — P-Xi \ as a basis for the calculation. 



3.1 Results for the X-distribution 

To make meaningful comparisons for the behavior of the X-distribution at different lattice 
spacings and different lattice sizes, it is necessary to fix the fraction / of the points that 
are examined on each lattice. In fact, it would be more appropriate to always speak of 
Xj-distributions. For a given low-lying mode we order the lattice sites by the magnitude of 
density and consider the top fV points, where V is the lattice volume. It is assumed that 
the underlying gauge field is strongest in the regions so selected [1]. For most of the results 
discussed here we fixed / = 0.1. This is mostly motivated by the fact that we intend to 
relate our findings to the ILM. Diakonov and Petrov [7] use the packing fraction i in their 
theoretical arguments. On the other hand, if one naively calculates the packing fraction as 
/ = nVp, where n is the density in fm^ and Vp the volume of a four-dimensional sphere of 
radius p, then one obtains approximately ^ if ILM values are used. We thus view f = Jq 
as a reasonable compromise. The qualitative and even quantitative conchisions that we will 
make do not depend on the precise value from the range quoted above as long as it is fixed. 

Our results for the X-distribution from the lowest two nonzero pairs of modes included 
for each configuration arc shown in Fig. 1. The histograms arc normalized so that the sum of 
the values in all bins adds up to unity. There are visible peaked maxima at X ±0.65, but 
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Figure 3: Comparison of X-distributions at /? = 6.55 for near-zero modes (solid line) and 
exact zero modes (dashed line) using the overlap Dirac operator. 



there appears to be very little change across the wide range of lattice spacings studied. To 
quantify this observation, we have calculated the average value of |X| with the histograms 
serving as a probability profile, i.e. 

< |X| > = J2px. (4) 

i 

where px, is the value at the bin with center Xj. For truly peaked distribution, this should 
reflect the approximate position of the peak. The results are shown in Fig. 2 indicating a very 
flat behavior as a function of the lattice spacing. The underlying dynamics thus does not 
appear to generate more chiral peaking closer to the continuum limit, nor do the positions 
of the peaks move appreciably closer to ±1. 

As we have emphasized in Sec. 2, the qualitative observation that the X-distribution 
exhibits peaked maxima at nonzero X is not sufficient for concluding that the instanton 
picture of topological charge fluctuations in the QCD vacuum is correct. However, the 
study of quantitative characteristics of the distribution might still be very illuminating. We 
illustrate this by comparing in Fig. 3 the X-distributions for the near-zero modes and for 
the exact zeromodes. The latter is the same distribution that would be expected in the limit 
of a very dilute instanton gas (or a dilute gas of arbitrary topologically nontrivial lumps) . It 
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Figure 4: Density of structures (in fm ^) as a function of the lattice spacing. 



is clear that, even at our smallest lattice spacing, the X-distribution of the near-zero modes 
is far less chirally peaked than that of the exact zero modes®, confirming the prevalent view 
that the dilute instanton gas picture is not realistic [7] . This also suggests that one rehable 
way to verify whether lattice QCD ^-distributions are consistent with the ILM prediction 
is a direct comparison to distributions from existing ILM ensembles. 

3.2 Structures in the Gauge Field 

Our next goal is to verify whether the observed X-distribution can be attributed to the 
underlying local structures with properties similar to those assumed in the ILM. If the 
low-lying mode arises from mixing of t'Hooft "would be" zero modes associated with ILM 
instantons, then this mode inherits the lumpy structure of the underlying gauge field, with 
lumps being of a prescribed size, shape and abundance. 

We identify such possible structures by finding the local maxima of density, d{n), in the 
mode, and requiring that the profile of density around them resembles a four-dimensional 
peak at least in some average sense. To be more precise, we start by locating the maxima 

^This is equally well reflected in the fact that we found typical values < |X| 0.52 (see Fig. 2), while 
at the same time < \X\ >= 1 for distribution strictly peaked at X = ±1. 
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over the distance -\/3, i.e. finding the set 

A4 = {n : d{n) > d{m) , \n — m\ < \/3 } (5) 
For these candidates, specified by their position n, we then compute the functions 

dn^^{r) = < d{m) > \„-m\=r n = ±1, . . . , ±4, r > (6) 

(n— m)./x>0 

representing the average of d{m) over the spherical shell of radius r centered at n, restricted to 
the points for which (n — m) has a component in the /^-direction. Local maximum is retained 
only if dn^^{r) decays monotonically from origin over the distance for all directions ^. 

Note that the motivation for choosing as a reference lattice distance in the above 
procedure is that it is a minimal distance for which the accidental occurrence of a structure 
can be statistically excluded for the largest lattice volumes we are working with. In other 
words, if we were to generate the density on sites of a 32^ lattice using random numbers, no 
structures would typically be found using the above prescription. Conversely, if a structure 
is identified in a low-lying mode, we take it as a signal that a nontrivial fiuctuation of the 
gauge field is causing its occurrence. 

Finally, from the structures identified in this way, only those were selected that have a 
chance of contributing to the peaks of the X-distribution. To make this final cut, we only 
retain structures for which the corresponding center n is among the 10% of points with 
highest density (/ = 0.1 was used to generate X-distributions), and for which the chiral 
angle satisfies |X(n)| > 0.5 as suggested by the position of the maxima in Fig. 1. The 
resulting density of structures as a function of lattice spacing is shown in Fig. 4, indicating 
a large disagreement with the ILM assumption at smaller lattice spacings. We will return to 
the interpretation of this behavior later. 

3.3 Sizes from Coherent Regions of Local Chirality 

To study the typical sizes of the structures, we first use the definition that is motivated 
only by the assumed lumpy behavior of the gauge field rather than the specific profile of an 
instanton. The mixing of "would be" zero modes in the lumpy background would produce 
regions of coherent local chirality in the near- zero modes, concentrated around the maximum 
of the lump. We thus define the radius of the structure at point n as the radius of the largest 
hypersphere centered at n, containing points with the same sign of local chirality, namely 

Rn = max { r : c(n) c(m) > 0, \n — m\ < r } (7) 

To assess the feasibility of this definition, we compute functions 

dnir) =< d{m) >|n-m|=r C„(r) =< C(m) >|n-m|=r (8) 

representing the average density and chirality over the spherical shell of radius r centered 
at n. In Fig. 5 we show the sample behavior of functions (i„(r)/(i„(0) and c„(r)/c„(0) for 
typical structures, and how this relates to the radius assigned by our definition (vertical line 
in Fig. 5). The examples of structures with high, medium and low intensity (density) are 
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Figure 5: The sample behavior of (i„(r) / (i„(0) (upper curve) and c„(r) / c„(0) (lower curve) for 
typical structures &t P — 6.55. The radius, of the corresponding region of coherent local 
chirality is shown as a vertical line. The figures are ordered so that the intensity of the peaks 
decreases vertically and the quality (degree of isolation) of peaks decreases horizontally. 
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Figure 6: Average radius, < i?„ >, of structures from regions of coherent local chirality. The 
lowest nonzero mode was used for calculation. Data for the three smallest lattice spacings 
were used to obtain the fit. The horizontal line represents the radius of an ILM instanton. 



put in the top, middle and bottom row respectively. Note that the level of "background" 
increases for lower intensities. Also, the structures with cleaner peaked behavior (typically 
more isolated ones) are put on the left and quality is decreasing to the right. The examples 
on the right represent the "worst" cases we have identified. From these examples it can be 
seen that the peaks of density and chirality are well contained within the determined radius. 

Using the definition (7), we have studied the average sizes of structures as a function of the 
lattice spacing with results shown in Fig. 6. The straight line is a fit using the three smallest 
lattice spacings considered. The average radius at a ~ 0.12 fm is close to the ILM value 
and consistent with a recent estimate at a similar lattice spacing [12]. However, it would be 
erroneous to conclude agreement with the ILM based on the data at single lattice spacing; 
the average value decreases significantly for finer lattices and the continuum-extrapolated 
estimate is in striking disagreement with the ILM. There appears to be a positive curvature 
in our data and a well-defined finite value (with rough estimate of about 0.15 fm) in the 
continuum limit. The fact that our procedure leads to a finite size in physical units is 
significant because it characterizes the physical size of regions of coherent local chirality that 
are necessary to develop the large rj' hairpin correlator required to solve the U{1) problem [1]. 
We also emphasize that while our definition of size was motivated by assuming the lumpy 



16 



Size distribution of Structures beta=6.55 V=32'^4 



J 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

Radius[fm] 

Figure 7: The distribution of sizes i?„ at /3 = 6.55. 

structure of topological charge fluctuations, it makes very good sense even if the structures 
do not carry the quantized values of topological charge, and are not related to instantons. 
The distribution of sizes a,t P — 6.55 is shown in Fig. 7. 

3.4 Sizes from Instanton- Assumed Profiles of Density 

If the near-zero mode of the Dirac operator is a mixture of "would be" zero modes associated 
with instantons, then peaks of the wave function should resemble the profile of a 't Hooft 
zeromode. In particular, in the ideal case we should have [3] 




where p„ is the radius of the instanton located at n. Thus the function, r^y„/(l — yn) 
is a constant (p^J independent of r for an instanton profile, and should be approximately 
constant for our structures if they represent the response to an instanton-likc fluctuation. 
However, we find that this is not the case and the shapes of the peaks in our low eigenmodes 
do not resemble the instanton ansatz. To illustrate this, we display the situation for a 
typical structure in Fig. 8. Shown are also fits to the instanton profile over the the distances 
0.00 — 0.06 fm, 0.06 — 0.12 fm and 0.12 — 0.18 fm. These fits are inconsistent with both the 
profile of the peak and with each other. Note that the shapes of the peaks are already 
averaged over all the directions, and should exhibit a robust behavior. 

17 



o 
O 



200 



150 



100 



50 



beta=6.55 config 1 struct 8 



beta=6.55 config 1 struct 8 



50 




5 10 15 20 25 30 35 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



r'^2 [lattice units] r[fm] 
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should be constant for an instanton profile. Right: The profile of the structure and attempted 
fits in the region 0.00 — 0.06 fm (leftmost curve), 0.06 — 0.12 fm (middle curve), and 0.12 — 
0.18 fm (rightmost curve). 



To see the inconsistency with the instanton profile on average, we choose a reference 
point Tref (in lattice units) and assign a radius to each structure through 



Pn(j'ref) — "^ref 



Unyref) ^^^-^ 



^ 1 - yn{rref) 

All of our structures are guaranteed to have a peaked behavior over the distance r^^j = 3 
in lattice units, but the vast majority of them decays over much larger lattice distances. 
To treat all the structures on the same footing we compute average radii from the above 
prescription for r'^^^ — 1, 2, . . . , 5 and plot the resulting dependence on lattice spacing in 
Fig. 9. Wide vertical bands at every lattice spacing reflect the fact that our structures 
can not be flt reasonably by the t'Hooft proflle. Nevertheless, the tendency towards sizes 
substantially smaller than the ILM value in the continuum limit is still obvious, thus showing 
that the overall conclusion from Fig. 6 remains valid even for this definition of the radius. 

We emphasize that contrary to results in Fig. 6, we do not assign much physical signifi- 
cance to the precise behavior of data in Fig. 9. The main purpose of showing these results 
is to reveal the marked inconsistencies with the instanton picture of the vacuum. 

3.5 Correlation Functions of Local Chirality 

We have also computed the chirality-chirality correlation function in the lowest Dirac nonzero 
mode V', namely 

Ccc(r) = {c{n)c{m))\n-m\=Ty'^ (H) 
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Figure 9: Average radius, < p„ >, of structures determined from the assumed instanton 
shape of the peak at different reference distances. Lowest nonzero mode was used for calcu- 
lation. 



This correlator has been studied in Ref. [12]. While the average information encoded in 
the correlator may be useful, it will not provide us with a detailed view of the space-time 
distribution of c(n). The point is that if this Ccdr^ has characteristic size p, one still cannot 
infer whether the typical space-time distribution of cin) in a configuration predominantly 
appears localized on spherical structures of radius p, or if it comes in one of the infinitely 
many other forms leading to the same average correlation. 

The average Cccir^jCcJ^^ from the lowest nonzero modes for our four ensembles are 
shown in Fig. 10. While larger statistics would be desirable, wc have estimated the "size" 
of the correlators by evaluating < r > and < >^/^ with C^At) used as a probability 
distribution. In every case we have cut off the integration at the distance where the correlator 
first turns negative. The resulting values are collected in Table 2. Note that values for < 
>^/^ are systematically higher than sizes obtained from coherent local chirahty (c.f. Fig. 6) 
which is to be expected for an inhomogeneous space-time distribution of fluctuations with 
various shapes. Nevertheless, the sizes of correlators we have obtained are still systematically 
lower than the ILM value for instanton radius. 
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Table 2: Estimates of the size for the correlator Ccc(^)- 
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3.6 The Interpretation of the Results 

The data in Figures 4 and 6 represent rather interesting new results and we would now like to 
elaborate on their interpretation. We have started from the assumption that the vacuum has 
local properties identifiable with ILM scenario, i.e. that it typically fiuctuates in such a way 
as to form (anti)self-dual lumps of approximately quantized topological charge. However, 
the identification of such presumed lumps resulted in quantitative characteristics that are 
in marked disagreement with the ILM on fine lattices. We observe many more structures of 
much smaller physical size (although this size remains finite even in the continuum limit). 
Given that the parameters of the ILM are rather tightly fixed [7], the true microscopic picture 
as seen by the lattice fermion appears to be very different from that envisioned by the ILM. 

Particularly puzzling from the instanton perspective is a large number of peaks that we 
observe in the lowest near-zero modes, implying that a propagating light quark feels many 
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more "kicks" from the regions of strong gauge field than one would expect based on the 
ILM. Indeed, assuming that the underlying gauge structures represent elementary instanton 
tunneling events between classical vacua leads to clear contradictions. For example from our 
results dX (3 = 6.55 one would estimate that the gluon condensate should be of the order 

< OIG^IO >^ 327r2n ^ 3271^ 50 fm"^ ^ 25 GeV^ where ?>2ti'^ is a contribution of a single 
instanton ^. This is to be compared with the accepted value < 0|G^|0 >~ 0.5 — 1.0 GeV^. 
Similar inconsistency arises from considering topological susceptibility under even weaker 
assumptions: Assuming that gauge structures are general lumps of quantized and approx- 
imately uncorrclatcd topological charge leads to the estimate of topological susceptibility 

< Q"^ > /V ^ n ^ 50 fm"'^ ^ 50(200MeV)*^, i.e. about fifty times the accepted value. At the 
same time, the topological charges of our lattices (determined as a byproduct of the overlap 
calculation) are in very good agreement with the expected value of topological susceptibility 
in the pure gauge vacuum ^. 

The above considerations strongly suggest that the topological charge contained within 
the structure is probably not quantized. One is led to the same conclusion from consideration 

''The estimates of this type are usuahy used to justify the ILM parameters [7]. 

®For example, from our small ensemble at /? = 6.55 we estimate the topological susceptibility x = 
0.92 ±0.50 fm"*. 



21 



of radii of regions of coherent local chirality. Our data indicate the value <itVi>Ri0.15fmin 

the continuum limit. However, the 't Hooft instanton formula for density of quantized charges 
valid at short distances, i.e. n{p) oc p®, is frequently invoked as a basis for concluding that 
the occurrence of instantons of this size should be essentially negligible (certainly relative to 
p 0.3 fm instantons for which the formula is still expected to be valid) [28, 7]. We interpret 
this marked mismatch as a manifestation of the fact that the underlying assumption for the 
't Hooft formula, namely the local qiiantization of topological charge, is probably not valid. 
Indeed, the suppression of small instantons is a consequence of their quantized topological 
charge: it is difficult to squeeze an entire unit of topological charge inside a small radius. 
However, without quantization, there is no reason to expect a suppression of smaller-scale 
fluctuations, which would simply contain proportionally less topological charge. 

Finally, we would like to return to Fig. 4 and discuss the rapidly increasing density of the 
structures with decreasing lattice spacing. While the density aX (3 = 6.55 is very large, the 
fraction of sites contained in the structures is still relatively small as shown by the packing 
fractions plotted in Fig. 11. The peaks thus retain their identity as can be seen from Fig. 5. 
However, let us elaborate for the moment on the hypothetical possibility that the density 
increases indefinitely, because this might suggest something unphysical even if topological 
charge is not locally quantized as we argued above. Indeed, how should one interpret a 
diverging density of structures in the continuum limit, especially if their average size (c.f. 
Fig. 6) remains finite? The point is that contrary to quantized topological charges, where 
the finite topological susceptibility forces the density to assume a finite physical value in the 
continuum limit, the peaks of non-quantized topological charge might well lose their identity 
closer to the continuum limit and their density be indeed unphysical. (However, the regions 
of coherent local chirality identified with the help of these peaks, and their sizes are still 
physical.) The picture that we have in mind is that of relatively isolated mountain peaks 
belonging to a larger mountain range and coalescing as the continuum limit is approached. 
Our algorithm to measure the size of regions of coherent local chirality would then sample 
the size of the mountain range rather than individual peaks. A more dynamical analogy is 
to imagine disturbing a calm surface of water. The regions where the water is above the 
original level (positive topological charge density) have various shapes and volumes, but on 
average have some typical size. Yet, within such regions there could be many local maxima 
at all length scales. While the sizes of individual maxima do not define the scale and may 
not be physically meaningful, the average radius measured relative to these maxima probes 
the typical finite size of these coherent regions. It is easy to imagine that in the continuum, 
the fiuctuations of topological charge resemble violently disturbed surface of water. In such 
a hypothetical scenario, the diverging density of structures would actually be very natural. 

We emphasize that in this paper we do not attempt to put forward any specific low-energy 
scenario for the behavior of topological charge fluctuations in the continuum limit. The 
purpose of the discussion in the previous paragraph is to argue that as soon as one abandons 
the local quantization of topological charge, then the data in Fig. 4 do not represent anything 
unexpected. 
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3.7 The Question of Dislocations 



In the lattice discussions of gauge field topology one has to face the possible problem of 
"dislocations" . While this notion is frequently used as an unspecified synonym for "lattice 
artifact", in the original discussion of Ref. [29, 30] it has a rather well defined meaning 
which wc will adopt. In particular, for a given lattice gauge action and given topological 
charge operator, it is the local structure in the gauge field with action smaller than the 
continuum action of an instanton, and with unit topological charge ^. In the loose sense, 
this is frequently pictured as a small instanton that is almost falling through the lattice 
and lives on 1-2 lattice spacings. Invoking entropy arguments, it was suggested in Ref. [30] 
that for the action-topological charge combination where this happens, it is possible that 
dislocations dominate the topological charge fiuctuations, and can lead to unphysically large 
(possibly diverging) susceptibility. This raises the question of whether the large density of 
structures that we see is caused by unphysical dislocations 

This question can be addressed at several levels. First of all, the arguments of Ref. [30] 
assume that the topological charge of the dislocation is concentrated on its core. Indeed, the 
idea that dislocations are important relies rather heavily on the picture of locally quantized 
topological charge. If the topological charge does not come in unit lumps, as we suggest 
here, then the issue simply does not arise. 

Secondly, it is important to realize that dislocation is a concept assigned to the pair 
gauge action-topological charge operator. While we work with Wilson gauge action which 
is supposedly susceptible to the possibility of dislocations, we use chiral fermions to measure 
topological charge. It has frequently been noted (see e.g. Ref. [7]), that fermions should not 
be sensitive to dislocations. This is indeed very plausible. As we have argued in Sec. 2.3, one 
naturally expects that infrared fermionic modes will be quite smooth and will not inherit 
the singular behavior of the underlying gauge fields. 

Finally, there is no hint of behavior symptomatic of dislocations in our data. For example, 
even the smallest of our structures are extended objects in physical units and span at least 4 
lattice spacings. Moreover, the estimates of topological susceptibility from topological charge 
measured by the overlap Dirac operator are in very good agreement with the accepted value 
of about 1 fm~^ (see footnote 8). If the underlying structures were dislocations (and hence 
carried topological charge close to unity), then this value should be substantially larger. 
This is one of the arguments that actually suggests that topological charge is not locally 
quantized, as we stressed in Sec. 3.6. 

Nevertheless, to eliminate any concern that dislocations may be a problem, we have 
performed a comparative study with a renormalization group improved gauge action. In this 
case, the possibility of dislocation dominance should be substantially reduced. In particular, 
we have employed the Iwasaki action [31], and studied overlap eigenmodes at two different 
lattice spacings. In the standard normalization for Iwasaki action, (see Refs. [31, 32]), we 
considered a 10^x12 lattice at f3 = 2.384 and a 14^x16 lattice at /3 = 2.592. This corresponds 
to lattice spacings a = 0.14 fm and a = 0.10 fm respectively, as determined from string 
tension by using the data of Ref. [32] and the two-loop beta function fit. 

The summary of our results is shown in Fig. 12. The structure of the resulting X- 

^We thank Tamas Kovacs for clarifying this. 
^°We thank Thomas Schafer for pointing this out. 
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Figure 12: Upper: X-distributions for two Iwasaki gauge ensembles considered. Lower 
left: Density of structures for Iwasaki and Wilson actions. Lower right: Average size from 
coherent regions of local chirality for Iwasaki and Wilson actions. 
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distributions is very similar to the structure for Wilson gauge ensembles as can be seen by 
comparing to Fig. 1. We thus applied identical procedures (with identical cuts) for finding 
the structures and determining their size as we did for Wilson gauge ensembles. Although the 
behavior of density for Wilson gauge action (see Fig. 4) might have suggested the possibility 
of dislocations, this is not the case because a similar pattern appears to be followed in the 
case of Iwasaki action as well. Similarly, the sizes of structures are systematically lower than 
the ILM values. It would obviously be interesting to investigate a larger window of lattice 
spacings with Iwasaki action, to be able to study the issues of continuum extrapolation. 
However, the purpose of this comparative study was to demonstrate that our conclusions 
are not sensitive to the choice of gauge action and that there is no problem of dislocations. 
Our results clearly indicate that this is indeed the case. 

4 Conclusions 

Uncovering the local structure of topological charge fluctuations promises to have profound 
implications for our understanding of low energy QCD. This expectation is born out of the 
fact that possible microscopic explanations for important phenomena such as spontaneous 
chiral symmetry breaking, the resolution of the Ua(1) problem, and the ©-dependence of 
QCD, are based on the picture of the vacuum wherein self-dual lumps of locally quantized 
topological charge (instantons) play a major role. It is consequently quite essential to deter- 
mine whether such picture (e.g. ILM) is indeed fundamental, or if there is in fact a different 
microscopic mechanism driving these important effects. 

Following upon the ideas of Ref. [1], we have explained here in detail why studying 
low-lying Dirac eigenmodes provides a natural and reliable approach for exploring the local 
nature of topological charge fluctuations (at least as implied by the ILM). The study of 
the X-distribution in low-lying modes has been designed as a tool for probing the local 
vacuum structure indirectly. While the available results [13, 14] on qualitative behavior of 
X-distribution could in principle be viewed as confirming the consistency of the instanton 
picture, we have pointed out here that the association with instantons is not unique. In other 
words, the double-peaked qualitative structure of X-distribution is a necessary condition for 
instanton dominance, but not sufficient. 

As a step towards resolving these issues, we have computed and analyzed low-lying modes 
of the overlap Dirac operator in Wilson gauge backgrounds over a wide range of lattice 
spacings (a « 0.04 — 0.12 fm). The double-peaked structure of the X-distribution has been 
observed with maxima at X ^ ±0.65 (using a fraction / = 0.1 of lattice points). The 
position of maxima shows very little change with lattice spacing without visible tendency to 
move toward the extremal values oi X = ± 1 , even though an exactly chiral lattice fermion 
is used. Our main initial aim in this study was to verify whether the characteristics of local 
structures contributing to the peaks of the X-distribution are consistent with quantitative 
characteristics of the gauge features in the backgrounds proposed in the ILM. In other words, 
we attempted to verify whether low-lying modes can (at least approximately) be viewed 
as mixtures of t'Hooft "would be" zeromodes generated by gauge structures resembling 
approximately p = i fm instantons of density n — 1 fm~^. To the contrary, we have found 
local characteristics indicating that the true nature of topological charge fiuctuations is very 
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Figure 13: The X-distributions for first 20 near-zero modes of configuration 8 at /5 = 6.2. 
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different from that envisioned in the ILM. 

We emphasize in this context that the underlying issue here is to determine whether 
one should view the ILM picture as a fundamental one, or as an effective description with 
phcnomcnological meaning only. For example, claims that the bulk of spontaneous chiral 
symmetry breaking is due to mixing of "would-be" zeromodes associated with instantons 
have only well defined microscopic meaning if the corresponding low-lying modes have the 
local structure predicted by this microscopic mechanism. To verify this proposition, we have 
no other choice but to examine these modes. This appears to be a well-defined problem 
with a well-defined answer. Our results indicate that the ILM scenario does not provide for 
accurate microscopic description of these modes and thus remains at the phenomenological 
level. It is this important distinction that we wish to stress in this paper. 

While the detailed arguments are summarized in Sec. 3.6, our main conclusion can be 
verified in an independent manner. To make this point qualitatively, we have calculated 
twenty near-zeromodes for configuration 8 from our j3 = 6.2 ensemble This configuration 
has (5 = 0, and the physical volume is such that according to the ILM it should contain 
about 3 or 4 instantons and antiinstantons. Hence, the subspace spanned by t'Hooft would-be 
zeromodes should have a dimension of that order. Consequently, it is the ILM prediction that 
chirally peaked X-distributions should be observed for 3 or 4 near-zeromodes (and possibly 
a few more) but the rest of the modes should resemble approximate free-field behavior with 
local chirality peaked around the origin. In Fig. 13 we plot the X-distribution at / = 0.1 
for all 10 pairs of near-zeromodes computed (the histogram is the same for both modes in 
a pair). Inspecting these results reveals that there arc at least 14 modes with significant 
double-peaked structure and none of them is peaked at the origin. As a matter of fact, at 
/ = 0.02 all 20 calculated modes exhibit the double-peaked behavior and this most likely 
persists even for higher modes. While this fact can hardly be explained by the ILM, it is 
not very surprising in view of the results presented here. On this particular configuration, 
we have identified 32 local structures. 

Even though we have concentrated on the ILM in this paper, we believe that our re- 
sults suggest a more general conclusion. In particular, as we have argued in Sec. 3. 6, it is 
very difficult (if not impossible) to reconcile the topological susceptibility of the pure gauge 
vacuum and the assumption that local structures in fermionic near-zeromodes are caused 
by underlying gauge excitations with locally quantized topological charge. This leads us to 
believe that the bulk of topological charge in the QCD vacuum is not locally quantized in 
integer units as suggested some time ago by Witten [6]. We will address this issue in detail 
in a forthcoming publication [33]. 

Acknowledgments: This work has been supported in part by U.S. Department of Energy 
under grants DE-FG05-84ER40154, DE-FG02-95ER40907, DE-FG02-97ER41027 and DE- 
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Kovacs. 

^^Needless to say, this was rather demanding on computer resources. 

^^Note that this observation makes also very unhkely the possibiUty that one could interpret several 
structures as forming a severely deformed quantized lump of topological charge which one could somehow 
still associate with the ILM instanton. In this case the subspace of would-be zeromodes should still be of 
size 3 or 4. 
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